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Abstract 

(n) 

We study properties of a non-Markovian random walk X} , I = 0,1,2, ... ,n, 
evolving in discrete time I on a one-dimensional lattice of integers, whose moves 
to the right or to the left are prescribed by the rise-and-descent sequences char- 
acterizing random permutations it of [n + 1] = {1,2,3, . . . ,n + 1}. We determine 

(n) 

exactly the probability of finding the end-point X n = Xn of the trajectory of such 
a permutation-generated random walk (PGRW) at site X, and show that in the limit 
n — > oo it converges to a normal distribution with a smaller, compared to the conven- 
tional Polya random walk, diffusion coefficient. We formulate, as well, an auxiliary 
stochastic process whose distribution is identic to the distribution of the intermediate 

(n) 

points X} , I < n, which enables us to obtain the probability measure of different 
excursions and to define the asymptotic distribution of the number of "turns" of the 
PGRW trajectories. 



PACS: 02.50.-r Probability theory, stochastic processes, and statistics; 
05. 40. -a Fluctuation phenomena, random processes, noise, and Brownian motion 



1 Introduction. 



Properties of random unrestricted or certain patterns-avoiding permutations of [n] = 
{1, 2,3, ... ,n} have been analyzed by mathematicians in group theory and combinatorics 
for many years [1-3]. Studies of several problems emerging within this context, such as, 
e.g., the celebrated Ulam's longest increasing subsequence problem (see, e.g., Refs.[4-6] and 
references therein), provided an entry to a rich and diverse circle of mathematical ideas 
[7], and were also found relevant to certain physical processes, including random surface 
growth [8-11] or 2D quantum gravity (see, e.g., Ref.[6]). 

In this paper we focus on random unrestricted 1 permutations from a bit different view- 
point addressing the following question: what sort of random walk one gets when random 
uniform permutations are used as a generator of the walk? Here we consider a simple 
model of such a permutation-generated random walk (PGRW), evolving in discrete time 
on a one-dimensional lattice of integers, in which model the moves of the walker to the 
right or to the left are prescribed by the rise-and-descent sequence characterizing each 
given permutation n = {tti, n 2 , 7T 3 , . . . , 7Tj, . . . , 7r n+1 }. In a standard notation, the "rises" 
(or the "descents") of the permutation n are such values of / for which iii < {iti > 
[12]. We note that such a generator is evidently different of those producing conventional 
RWs, since here a finite amount n + 1 of numbers is being shuffled and moreover, neither 
of any two numbers in each permutation may be equal to each other; this incurs, of course, 
correlations in the rise-and-descent sequences and implies that the resulting PGRW is a 
non-Markovian process. 

In this paper we determine exactly many characteristic properties of such a random 
walk. First, we define the probability V n (X) of finding the end-point X n of the PGRW 
trajectory at site X. We show next that in the long-time limit V n (X) converges to a 
normal distribution in which the effective diffusion coefficient D = 1/6 is three times 
smaller than the diffusion coefficient (D = 1/2) of the conventional Polya random walk 
in one dimension (ID), which signifies that correlations in the generator are marginally 

1 In the sense that all patterns are permitted. 
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important. Indeed, assuming that correlations in the sequence of rises depend only on 
the relative distance between their positions, we deduce from V n (X) the two- and four- 
point correlation functions explicitly and demonstrate that correlations extend effectively 
to nearest-neighbors only. Next, at a closer look on the intrinsic features of the PGRW 
trajectories, we formulate an auxiliary Markovian stochastic process. We show that, despite 
the fact that this process is Markovian, while the PGRW is not a Markovian process, the 
distribution of the auxiliary stochastic process appears to be identic to the distribution of 
the intermediate states of the PGRW trajectories. This enables us to obtain the probability 
measure of different excursions of the PGRW, deduce a general expression for /c-point 
correlation functions in the rise-and-descent sequences, as well as to evaluate the asymptotic 
form of the distribution of the number of "turns" of the PGRW trajectories, i.e. the number 
of times the walker changes the direction of its motion up to time n. Finally, we discuss, 
in the diffusion limit, the continuous space and time version of the PGRW. 

The paper is outlined as follows: In Section 2 we formulate the model more specifi- 
cally. In Section 3 we derive the distribution function of the end-points of the walker's 
trajectories exactly and analyse its asymptotical behavior. In Section 4 we consider the 
correlations in the rise-and-descent sequences and obtain explicit results for two- and four- 
point correlation functions. Further on, in Section 5 we introduce an auxiliary stochastic 
process which has the same distribution as the intermediate points of the walkers trajec- 
tory X^ n \ I — 1, 2, . . . , n — 1, and obtain the probability measure of different excursions. 
In Section 6 we determine the asymptotic form of the distribution function of the num- 
ber of "turns" of the PGRW trajectory, relating it to Appoint correlation functions in the 
rise-and-descent sequences. Next, in Section 7 we discuss the continuous space and time 
analogue of the PGRW. Finally, in Section 8 we conclude with a brief summary of our 
results and discussion. 
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2 Model. 



Let 7T = {/Ti, 7r 2 , 7T3, . . . , 7r ra+ i} denote a random, unconstrained permutation of [n + 1]. We 
rewrite it next in two-line notation as 



(1) 



1 2 3 ... n + 1 

7Ti 7r 2 7T 3 ... 7T n+ i 

where the numbers in the first line will be thought of as the values attained by a discrete 
"time" variable /. We call, in a standard notation, as the "rise" (or the "descent") of the 
permutation n, such values of I for which i\i < 7r i+1 {jii > 7r i+1 ) [12]. 

Consider now a walk evolving in time / on an infinite in both directions one-dimensional 
lattice of integers according to the following rules: 

-at time moment I = the walker is placed at the origin. 

-at time moment I = 1 the walker is moved one step to the right if 7i"i < 7r 2 , i.e., I = 1 
is a rise, or to the left if i\\ > 7r 2 , i.e / = 1 is a descent. 

-at time moment / = 2, the walker is moved one step to the right (left) if 7r 2 < 7r 3 
(7r 2 > 7r 3 , resp.) and etc. 

Repeated / times, this procedure results in a random, permutation-dependent trajectory 
X\ n \ (/ = l,2,...,n), 

*, (B) =I>, ( 2 ) 



k=i 



where the "spin" variable is 



Sk = 6{^k+l ~ 7T fc ), (3) 

6{x) being the theta- function with the properties 

{+1 if x > 0, 
(4) 
—1 otherwise. 

The questions, which we address here, are a) the form of the probability V n (X) of finding 
such a walker at position X at time moment n, b) asymptotical behavior of V n (X), c) 
the probability measure of different excursions X^ and d) the distribution function of 



the number of "turns" of the PGRW trajectories. As a by-product of our analysis, we 
also determine some specific correlations in permutations, embodied in the moments of 
the walker's trajectories., i.e. the fc-point correlation functions in the rise-and-descent 
sequences, as well as the asymptotic distribution of the number of peaks and throughs in 
random permutations. 

3 The probability distribution. 

One possible way to determine V n (X) is to reconstruct it from the moments < X^ q >, 
(q = 1, 2, 3, ... ), of the end-point X n of the walker's trajectory, X n = X^ n = Ylk=i Sfc ' 



where the angle brackets denote "averaging" procedure - summation with respect to the 
set {it} of all possible permutations n, weighted by their total number (n + 1)!. 

Pursuing this kind of approach, we may then represent the summation over {n} as a 
multiple summation over states of a set of independent variables {a«}, each running from 
1 to n + 1: 



where the sum with the subscript j) extends over all i, j pairs such that i, j G 1,2,3,..., n+ 
1, excluding i = j, 5 ai>aj is the Kroneker- symbol, while the factor 



accounts for the constraint that 7Tj ^ ttj for any i and j, i ^ j. 

Therefore, calculation of moments of X n amounts to the evaluation of (rather) special 
correlations in an (n + l)-state Potts-like model with long-range interactions on the ID 
lattice containing n + 1 sites, which may be approached, in particular, using a non-trivial 
transformation of variables proposed in Ref. [13]. It appears, however, that within this 









4 



approach the calculation of already the second moment < X\ > is quite laborious, although 
the final result is simple. 

On the other hand, V n {X) can be obtained straightforwardly, in an explicit form, 
following a different way of reasonings. To do this, let us recall that according to the 
definition of the PGRW, at the Z-th step the walker makes a move to the right (left) if / 
is the rise (descent) of the random permutation n. Now, let A/j (A/j,) be the number of 
"rises" ("descents") in a given random permutation ir. Then, for this given permutation 
7r the end-point X n of the walker's trajectory is just 

X n =Afi-AT i , (8) 

which can be rewritten, using an evident " conservation law" A/j- + A/j = n as 

X n = 2A/j - n. (9) 

Therefore, for a given permutation n, the end-point X n of the walker's trajectory xj n ^ is 
entirely fixed by the value of the number of rises in this permutation. 

Now, a total number of permutations n of [n + 1] having exactly A/j rises is given by 
the so-called Eulerian number [14]: 

/ \ A/t+l 

f;H ( - ir (") w+i -' r 

where (°) denotes the binomial coefficient. Consequently, in virtue of Eq.(9), we have that 
the probability V n (X) of finding the walker at position X at time n, (—n < X < n), is 
given by 

X 2 1 

Using the representation of the Eulerian number in Eq.(lO), the result in Eq.(ll) can be 
also written down in explicit form as the following finite series: 

^>^THr(»r)(^- r )"' m 



We derive next several useful integral representations of the distribution function V n {X) 
and of the corresponding lattice Green function [15], which will allow for an easy access to 
the large-n asymptotical behavior of V n (X). 

Making use of the equality [16]: 

[(n+2+6)/2] . . r,\ / . o . 7 \ n+1 poo / ■ n \\ n + 2 



r=0 



(13) 

we find that V n {X) admits the following compact form: 

V .(X) ji±i^mj- (?my +2 cos(xk)dk . (14) 

Note that the integral representation in Eq.(14) is different from the usually encountered 
forms since here the upper terminal of integration is infinity. As a matter of fact, an 
integral representation of V n {X) in Eq.(ll), in which the integration extends over the first 
Brillouin zone only, has a completely different form, compared to that in Eq.(14), and 
reads 

= jdw, 1 r" + {k) s^ com ) cos(xk)dk = 

= *»>(n+l)l £ " + (?)] S " 1 " +2( ' ;) C0S < A '*)*' < 15 > 

where ^ n (k) is the polygamma function [16]. Derivation of the result in Eq.(15) is outlined 
in the Appendix A. 

Consider finally the asymptotic forms of the probability distribution function. To do 
this, it is expedient to turn to the lattice Green function Q(X,z) associated with the 
distribution function of the end-point of the PGRW trajectories. From Eqs.(14) and (15) 
we find that the lattice Green function of the end-point of the PGRW is given by 



,s m 2 (k) 



1_ r sin (z sin(fc)) cos(XA:) 
xz J : :! ■ 



nz J Q sin (k — zsm(k)) 
6 



Now, in the limit z — > 1 , (which corresponds to the large-n behavior of the distribution 
function V n (X)), one finds that the leading behavior of G(X, z) is as follows: 

g(x, * -> i-) ~ ^Y^y exp {-VW^lx]) , (17) 

which implies that in the asymptotic limit n — > oo, the probability V n (X) of finding the 
end-point of the walker's trajectory at position X converges to a normal distribution: 

^M^F^HB- (i8) 

Therefore, in the limit n — > oo the correlations in the generator of the walk - random 
permutations of [n + 1], appear to be marginally important; that is, they do not result 
in anomalous diffusion, but merely affect the "diffusion coefficient" making it three times 
smaller than the diffusion coefficient of the standard ID Polya walk. We will consider the 
form of correlations in the next section. 



4 Correlations in the rise-and-descent sequence. 

We address now a somehow "inverse" problem: that is, given the distribution V n (X), we 
aim to determine two- and four-point correlations in the rise-and-descent sequences. 

We start with the analysis of two-point correlations. From Eq.(5) we represent the 
second moment of the walker's displacement as 



<^) = <(A/- T -A/i) 2 ) 



n 

.1=1 



n—1 



= n 



(19) 



2 E E 

h=i j2=ji+i 

where we denote as Cj^ J2 the two-point correlation function of the rise-and-descent se- 
quence, 

41 = (*»*»)■ (2°) 
From our previous analysis, we know already that (X%) ~ n/3 as n — > oo, (see Eq.(18)). 
Consequently, in this limit the sum on the rhs of Eq.(19) is expected to behave as 



n-l 



E E c 

31=1 J2=jl + 1 



(2) 
'31 ,32 



n 
3' 



(21) 



where the sign "-" signifies that the rise- and- descent sequence is predominantly anticorre- 
lated and the probability of having two neighboring rises (or two descents) is lower than 
the probability of having a rise neighboring to a descent. The question now is how do the 
correlations decay with a relative distance m = \j 2 — ji|? We will answer this question here 
assuming that Cj^- 2 is function of the distance \j 2 — ji\ only, i.e. 

^ 2 -^= C(2) (-=b 2 -Jil)- (22) 

We note that such an assumption seems quite plausible at the first glance, since the rises 
(and descents) are evidently uniformly distributed on the interval [l,n]. We will show in 
the next section that it is actually the case using an auxiliary stochastic process having 
the same distribution of the intermediate steps as the PGRW. 
Introducing the generating function of the form 

oo 

*< 2 >(s)=X; (23) 

n=0 

we get from Eq.(19) the following relation: 

* ,2, m = (T^+(tV 2,w ' (24) 

where (z) is the generating function of two-point correlations: 

oo 

C^\z) = ^z m C i2 \m). (25) 

m=l 

Consequently, in order to define C^(m), we have to evaluate X^\z). To do this, we 
proceed as follows: define the Fourier-transformed V n (X) as 

n 

V n (k)= exp(ikX)V n (X). (26) 

X=-n 

Noticing first that V n {n + 1) = 0, then, making use of the representation in Eq.(ll) and 
changing the summation variable, we get 

^ )= ^zMf exp(2ttr) ^ i y (27) 
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Further on, using the equality 

n+l 



E 

r=0 



n+l— r 



n + l 



1 -y) n+2 Li_ n _i(y), 



(28) 



where Li_ ri _ 1 (j/) denotes the polylogarithm function [16], we obtain, by setting y = exp(— 2ik), 

(2i) n+2 sin n+2 (A;) 



™ (n + l)l 

Next, taking advantage of the expansion [16] 



Li_ n _W exp(-2ik) . 



(29) 



Li_ n _ 1 (y) = (n + 1)! In 



y 



—n—2 oo 



-E 

j"=o 
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n+2+j 



j\(n + 2+j) 



(30) 



where Bj stand for the Bernoulli numbers, we arrive at the following result: 



V n {k) 



f sin(/c) 



n+2 



with 



and 



(_l)(«+l)/2(2A;)»+3 

(n+l)! 



1 - 



ji(k) for n odd, 
j 2 (^) for n even, 



(31) 



E 

3=0 



[-l) j B. 



n+3+2j 



(2j + l)!(n + 3 + 2j) 



2.7 



-lf/ 2 (2A;) n + 2 g 



n+2+2j 



(32) 



(33) 



(™ + !) ! j^(2j)\(n + 2 + 2j) 

Now, we aim to determine an d explicitly Expanding V n {k) defined by 

Eqs.(31), (32) and (33) into the Taylor series in powers of k, we find that V n {k) obeys 

V n {k) = l-^(n + 2-25 nfi )k 2 + 

+ 3^ ((5n + 8) (n + 2) - 165 n , - 24<5 n>1 + 85 n , 2 ) A; 4 + O (k e ) , (34) 



which yields 



and 



X £) = l -(n + 2-25 nfl ) 



= ^ ((5n + 8) (n + 2) - 165 n , - 245 n>1 + 85 n , 2 ) . 



(35) 



(36) 



We note parenthetically that the obtained distribution V n (X) in Eq.(ll) appears to be 
platykurtic, since the kurtosis excess 



x n} 6 / 1 



j£2 \ ~ ,n V "~ 



is negative. 

Now, from Eq.(35) we find 



* (2) W = ^i (38) 
3(1-2;) 



which implies that C^ 2 \z) = —z/3 and hence, that 

I — if m = 1, 

C^\m) = l 3 (39) 

[0 ifm>2. 

Note that Eq.(39) signifies that for m > 2 the two-point correlations C^(m = \j 2 — in 
the rise-and-descent sequence decouple into the product < >< Sj 2 > and hence, vanish. 
Consequently, two-point correlations are short-ranged extending to the nearest neighbors 
only. Therefore, the probability p^f (m) of having two rises (or two descents) a distance m 
apart of each other is 

m (m)= m (m)^ 1 + C f (m) = <^ (JO) 





while 

puM = PiA m ) = = { f (4i) 

Consider next behavior of the four-point correlations: 

^31,32,33,34 = { t S h S j2 S j3 S 34^ (42) 

where ji < j'2 < J3 < ji- From our result in Eq.(39) it seems natural to expect that Cjfj 2 j 3 j 4 

vanishes as soon either (or both) j 2 — ji or j 4 — j 3 are greater than unity. Consequently, 
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only non-vanishing four-point correlations are of the form 

C (4) (m) = (s^Sj^S^+m^S^+rn^y (43) 

To evaluate C^(m), we will proceed along the same lines as we did with the analysis of 
two-point correlations. From Eq.(5) we have 

<X 4 > = <(M-Ar i ) 4 > = 



n 

.1=1 



jWl^ . . . g "V, 



E\ 1 2 
^-rn r> ( 44 ) 
m 1 !m 2 ! • • -m n ! 

where the sum extends over all positive integer solutions of equation mi+m 2 + - • •+m n = 4. 
Taking into account our results for C^(m), as well as noticing that sf — sf = 1, while 
sf = Si for any I, we find that 

< x «> = l n + 4>r^ e{n ~ 2) + 2 3TiT (n " ^dM" " 2 ) + 



n— 3 

- 3) + 2) C (2) (l)#(n - 3) + 4! ^ (n - 2 - m)c (4) (m)fl(n - 4). (45) 

m=l 

Multiplying both sides of the last equation by z n and performing summation, we find the 
following relation between the generating function X^(z) of the fourth moment of the 
walker's displacement and the generating function (z) of the four-point correlations: 

On the other hand, from Eq.(36) we have that 

(4) z(35z+15-80z 2 + 48z 3 -8z 4 ) 

x {z) = w^w ' m 

which yields, eventually, 

c<4, W = |i^) = ^ + 5( z2 + z3 + 24 + -)- < 48 > 

and hence, 

C^(m) = {}$ (49) 
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Note that Eq.(49) signifies that the four-point correlations decouple into the product of 
nearest-neighbor two-point correlations, C^(m) = C( 2 )(1)C( 2 )(1), for m > 2. 

Finally, noticing that three-point correlations of the form (^SjSj + iSj + 2 are equal to zero, 
we may straightforwardly calculate the probabilities of several particular configurations 
involving three and four rises and descents. In particular, the probability of having two 
neighboring rises and a descent at distance m apart of them is given by 

Ptt,l( m ) = ^(( 1 + Sj)( 1 + Sj+l)( 1 -Sj+m+l)) = 




if m — 1, 

8 (50) 
- ifm>2, 

while the probability of having two neighboring rises and another rise at distance m apart 
of them obeys: 



(51) 



In a similar fashion, we get that the probability of having two rises and another pair of 
rises (descents) at distance m apart of them is given by 

if m — 1, 

120 (52) 

36 ifmS2 ' 
the corresponding probabilities involving three rises obey 

{— if m — 1, f — if m = 1, 

18 ifm>2, {- ifm>2, 

while the probabilities involving two rises and two descents follow 

— if m — 1, f — if m = 1, 

} 5 PTT,u( m ) = \ 2 ]° PlT,uM = 

- if m > 2, — if m > 2, 
9 ^36 ~ 
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Equations (50), (51) to (54) show explicitly that the probabilities of different configurations 
of the rise-and-descent sequences depend not only on the number of rises, but also on their 
order within the sequence. 

5 Trajectories Xj for 1 < n. 

So far, we have defined the distribution of the end-points of the trajectories X^ n \ but, 
apart from the explicit results on the form of two- and four-point correlations, which allow 
us to reconstruct trajectories X[ n) with n = 4 (see Fig.l), we do not have an access to 
information on the distribution of the intermediate points I — 1, 2, 3, . . . , n — 1. On the 
other hand, it seems to be quite non-trivial. Indeed, we deal with a random walk, which 
makes a move of unit length at each moment of time with probability 1, but nonetheless 
looses somehow two thirds of the diffusion coefficient. Given that the correlations are short- 
ranged and extend to nearest-neighbors only, one might expect that the PGRW behaves 
effectively as an " antipersistent" RW [15] with a short-range one-step memory, such that 
the probability, in view of the results in Eq.(40) and (41), of making on the Zth step a 
move in the same direction at which it made a move on the (/ — l)st step is 1/3, while 
the probability of changing the direction is 2/3. One readily verifies that for such a walk 
X^j ~ n/2, i.e. the reduction in the "diffusion coefficient" is smaller than the one we 
actually find for the PGRW, ((^X^j ~ n/3). As a matter of fact, as we observe in Fig.l, 
for the PGRW the memory of the " antipersistency" is stronger and depends not only on 
the number of steps to the right or to the left, which the walker has already made, but also 
on their order. In other words, the PGRW represents a genuine non-Markovian process. 

5.1 The probability distribution of x[ n ^ for 1 < n. 

To determine the structure of excursions X\ n) of the PGRW we adapt a method proposed 
by Hammersley [17] in his analysis of the evolution of the longest increasing subsequence, 
and elucidated subsequently in Ref.[18]. The basic idea behind this approach, which we 
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exploit here, is to build recursively an auxiliary Markovian stochastic process Yi, which is 
distributed exactly as xj- n \ 

At each time step /, let us define a real valued random variable xi+i, uniformly dis- 
tributed in [0, 1]. Consider next a random walk on an infinite in both directions one- 
dimensional lattice of integers whose trajectory Yj is constructed according to the following 
step-by-step process: at each time moment I a point-like particle is created at position xi + ±. 
If xi + i > xi, a walker is moved one step to the right; otherwise, it is moved one step to the 
left. The trajectory Yj is then given by 

i 

Yi = ^20(x k+1 - x k ), (55) 

k=l 

where 6{x) is the theta-function defined in Eq.(4). We note that the joint process (a^+i, Y t ), 
and therefore Y h are Markovian since they depend only on (xj,Yj_i). Note also that Y\ 
is the sum of correlated random variables; hence, one has to be cautious when applying 
central limit theorems. A central limit theorem indeed holds for the Markovian process YJ, 
but the summation rule for the variance is not valid. 

Now, we aim to prove that the probability V{Yi = Y) that at time moment / the 
trajectory Y\ of the auxiliary Markovian process appears on the site Y has exactly the 
same form as the Eulerian distribution obtained for the end-point of the trajectory x[ n \ 
That is, V{Yi = Y) obeys 

\l + (-l)Y+i] /l + l\ 
X 2 1 

To prove this statement, we assign a random permutation n of [n + 1] to each realization 
of sequence {xi}, / = l,2,...,n+l, by ordering the x/-s and requiring that tti is the row 
of xv. 

X Wl < X W2 < .. < 2-7T n +l (^^) 

Next, we note that according to the definition in Eq.(55) each random trajectory Y/, / < n, 
is entirely determined by ordering of the corresponding sequence {xj}j</ + i. This implies 
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X 



Figure 1: A set of all possible PGRW trajectories for n — 4. Numbers above the 
solid arcs with arrows indicate the corresponding transition probabilities. Dashed-lines 
with arrows connect the trajectories for different values of I. 



that Y\ is adequately determined by the permutation it defined by such an ordering proce- 
dure. 

One writes next the weight p(Yi) of a given trajectory Y% as an integral over the se- 
quences {xi}i<i + i generating Y h which can be represented as a sum over all corresponding 
permutations: 

p(Yi) = / dxi . . . dx n+1 = ^2 Pi 71 )' ( 58 ) 

J{x t } generates Yi n generates Yl 

where p(tt) is the probability of a given random permutation ir. On the other hand, p(tt) 
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obeys: 

dxi . . . dx n+1 = J dx n+1 J dx n .. J dx 1 = - (59) 

X 7T - 1 [7T ( k ) - 1] — X k ^ X TT ~ 1 [TT (fc) + 1] 

where x^--i(o) = 0, x vr -i( n+2 ) = 1 and ff^ 1 denotes the functional inverse of ir, such that 
7r _1 (7Tfc) = fc. Finally, noticing that integration over the variables £j>/+i obviously gives 
1, in view of the purely iterative definition of the process YJ, which is independent of the 
"future" Xi>i+i, we have 

E PW= E J^hy- (60) 

7r generates Yj 7r generates Yj 

which is precisely the probability of the trajectory with I < n. 

Therefore, Y t and x\ n) are identically distributed. We emphasize that the distribution 
of Yi, and therefore the distribution of X^ n \ do not depend on n for I < n, which is a rather 
counter-intuitive result for the process xj n \ Consequently, the probability V{X^ = X) 
that at any intermediate step I, I — 1, 2, 3, . . . , n — 1, the trajectory x\ n ^ appears at the 
site X obeys 

V(x\ n) =X) = Vi(X), (61) 

where Vi(X) is defined by Eq.(ll). 

We note here parenthetically that the result in Eq.(61) allows us to study the properties 
of the random process X^ which depend on the intermediate states, such as, e.g., the 
span, the maximal excursions within time n, the time spent on positive sites and so on. 
In particular, it might be instructive to compare the growth of the average length L n 
of the longest subsequence of (not necessarily consecutive) rises, i.e. of the so called 
"longest increasing subsequence" , and the growth of the average maximal positive excursion 
< max{X n } > of the PGRW, which is also supported by rises in random permutations. 
Since we realized that the PGRW converges in distribution to standard Polya walk with 
diffusion coefficient D — 1/6, we may estimate < max{X n } > as [15]: 

1/2 / o \ V 2 



/4D\ 1 ( 2 \ 1 

< max{X n } >~( — j • y/n = f — J • y/n w 0.461 • ^/n 



(62) 
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Hence, the growth of the maximal positive displacement of the PGRW proceeds at a slower 
rate than that of L n , L n ~ 2y/n, due to a numerical factor, which is more than four times 
smaller. 



5.2 Probability measure of different trajectories. 

The equivalence of the processes Yi and xj 1 ^ allows us to determine explicitly the prob- 
ability of any given trajectory. We note that in the permutation language, this problem 
amounts to the calculation of the number of permutations with a given rise-and-descent 
sequence and has been already solved in terms of an elaborated combinatorial approach in 
Refs. [19-21]. Here we propose a novel solution of this problem. 

Let X\ n) be a given trajectory generated by a random uniform permutation n of [n + 1]. 
This trajectory, according to the definition of the PGRW, can be uniquely defined in terms 
of the sequence of rises (f) and descents (I) characterizing the permutation n. Now, let i\ 
and /j denote the integral operators of the form 



7| = / dx • and i\ = / 
Jx Jo 



dx ■ (63) 



Further on, to each /-step trajectory xj n \ we associate a polynomial Q x ( n )(x) G R[x] of 
degree /, defined as 

i 

Q x (n) (x) = Y[ii-l = hhh ■■■Ii- 1, (64) 
i=i 

where Jj assumes either of two values - 1^ and prescribed by the direction of the "arrow" 
at the i-th step in the corresponding sequence. In particular, the polynomial Q x (n) (x) 
corresponding to the / = 5-step trajectory X^ = {T, T, I, T, T} will be 



— hhhhh ' ^ — 








rx2 




= / dx\ I dx 2 




dx 3 


J X J X\ 


10 


J XZ 


O Ju Ju 


4 

X 


x 5 


40 8 12 


24 


^ 120' 



1 

dxd I dx 5 ■ 1 = 

X4 



(65) 
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Then, the desired probability measure of a given trajectory xj^ is given by 



p(X< B >) = f 
Jo 



Q x (n) (x)dx 



■ 1 



(66) 



i=i 



Note that this probability measure is not homogeneous, contrary to the measure of the 
standard Polya walk. 

The result in Eq. (66) may be compared with the analogous expression for the number of 
permutations with a given rise-and-descent sequence found by Niven [20]. Following Niven, 
consider a fixed up-and-down arrow sequence of length / and denote by k±, k 2 , ■ ■ ■ , k r the 
positions of downarrows (r is the total number of downarrows) along the sequence. Suppose 
next that this up-and-down sequence corresponds to some random permutation ix of [n+ 1] 
such that, according to conventional notation, an -up-arrow represents a rise, while a down- 
arrow is a descent. A question now is to calculate the number iV(X^) of permutations 
generating a given up-and-down sequence (or, in our language, a given trajectory X^). 
Elaborated combinatorial arguments show that N(X^) equals the determinant of a matrix 
A of order r + 1 whose elements a it j (where i stands for the row, while j - for the column) 
are binomial coefficients ( k hi J, where k = 0, k r+ i = 1 + 1, and (™) = if n > m [20]. 
Consequently, an alternative expression for the probability measure p{x[ n ' > ) may be written 
down as 



p{X\ 



(/ + !)! 



det 



(i 1 ) 



(fe) 
(&) 



,k 2 , 
l 



V 1 (V) ft 1 ) (V) ft 1 )/ 

One may readily verify that both Eq.(66) and Eq.(67) reproduce our earlier results in 
Eqs.(50) to (54) determining the probabilities of different four step trajectories. 

Finally, similarly to our Eqs.(20) and (42), we may define a general, /c-point correlation 



ft 
ft 



ikl) ' 

(I 2 ) 

(kr) 



(67) 



function , of the rise-and-descent sequence: 



Cju..,j k — { S ji S j'2 S j:i ' ' ' s j k / - 



(68) 
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where the "spin" variable Sk has been defined in Eq.(3). In terms of the auxiliary process 
Yi, Cj^ j can be rewritten formally as 

Ik \ „i pi k 

C hldk = ( II 9 ( x ii+ 1 ~ x u) ) = / dx h ■ ■ dx ik+i II 9 ( x ii+ 1 ~ x ii)i ( 69 ) 
W /{.*} J ° J ° 

where the brackets with the subscript {x^ denote averaging with respect to the distribution 

of the ensemble of variables Xk- 

We note that the purely nearest-neighbors nature of the fc-point correlations, (which we 

have taken for granted in Section 4), indeed becomes quite transparent for the process Y h 

such that Cj^ j k factorizes automatically into a product of the corresponding correlation 

functions of the consecutive subsequences, in which all jk differ by unity, as soon as any 

of the distances jk+i — jk exceeds unity. On the other hand, the correlation function 

Q(k) _ C^ k ) +1 - +k of a consecutive sequence of arbitrary order k can be obtained recursively 

using Eq.(69). To do this, we first note that C {k) can be represented as 

C {k) = fdxMx), (70) 
Jo 

where fk(x) are polynomials of order k, which obey the recursion 

fk(x) = - dxf k -i(x) + / dxfk-i(x), f (x) = I. 

J0 Jx 



(71) 



From this recursion, one finds immediately that are defined through: 

(-i) fc : 
(fc + i 



C W = ^LH!^ C ( fc - 1 -^, C(°) = l, (72) 



p=0 

such that the generating function of obeys: 



CW(z) = f2C W * = ^^-- (73) 



k=0 



Consequently, are related to the tangent numbers [22] , and are given explicitly by 

C< ' = ( t + 2)l fc+2 ' < 74 > 



19 



where are the Bernoulli numbers. Note that since the Bernoulli numbers equal zero for 
k odd, we find that correlation functions of odd order vanish, i.e. C( 2fc+1 ) = 0. 

Finally, using the asymptotic expansions for the factorial and the Bernoulli numbers, 
we get that for k >> 1, the /c-point correlation functions (k is even) decay as 

C W „ 2{-l) k ' 2 (=-) , (75) 

i.e. the /c-point correlations fall off exponentially with k with a characteristic length = 
ln- 1 ( 7 r/2). 



6 Distribution of the number of "turns" of the PGRW 
trajectories. 

In this section we study an important measure of how scrambled the PGRW trajectories 
are. This measure is the number M of "turns" of an n-step PGRW trajectory, i.e. the 
number of times when the walker changes the direction of its motion. We focus here on the 
asymptotic form of the distribution V(M , n) of the number of "turns" of an n-step PGRW 
trajectory. 

In the permutation language, each turn to the left (right), when the walker making 
a jump to the right (left) at time moment j jumps to the left (right) at the next time 
moment j + 1 corresponds evidently to a peak T! ( a through J,f) of a given permutation 
n, i.e. a sequence Kj < itj+i > itj+2 > < 7Tj +2 ). Consequently, the distribution 
function V(Af,n) of the number of "turns" of the PGRW trajectory, (i.e. the probability 
that an n-step PGRW trajectory has exactly J\f turns), is just the distribution of the sum 
of peaks and throughs. The latter can be defined apparently using the so-called peak 
numbers P{n + l,m) of Stembridge (i.e., the number of permutations of [n + 1] having m 
peaks), which obey the following three-term recurrence [23]: 

P(n + 1, m) = (2m + 2)P(n, m) + (n + l- 2m)P(n, m-l). (76) 
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Below we will evaluate the asymptotic form of such a distribution function using a different 
type of arguments, other than the Stembridge formula in Eq.(76), expressing V(Af,n) 
through the correlation functions of the rise-and-descent sequences. 

Let J\f p and Af t denote the number of peaks (i.e., the number of configurations iij < 
7ij + i > 7Tj +2 ) and throughs (i.e., the number of configurations iTj > tyj +1 < 7Tj +2 ) in a 
random permutation of [n + 1] , respectively. These realization-dependent numbers can be 
written down explicitly as 

^ n— 1 

and 

_^ n— 1 

M = ^(l-^)(l + ^i), (78) 

where the "spin" variable Sj has been defined in Eq.(3). The total number of turns Af of 
the RW trajectory, generated by a given permutation, is then given by 

n—l 

N = N p + N t = - (1 - Sj8 j+1 ) (79) 
j'=i 

We seek now the characteristic function Z n (k) of A/", defined as 

Z n (fc) = ( exp (ikAf) ) = exp ^ exp ^ ^ SjSj+ ^j ^ (80 ) 

where the angle brackets denote averaging with respect to random permutations of [n + 1]. 
Note that Z n {k) in Eq.(80) can be thought of as a partition function of a somewhat 
exotic one-dimensional Ising-type model in which the "spin" variables Sj of different sites 
are functionals of consecutive numbers in random permutations and thus possess rather 
specific correlations. 

Now, since SjSj + i = ±1, the characteristic function Z n {k) in Eq.(80) can be written 
down as 



Z n {k) = 




where t = tanh(ik/2). Averaging the product in Eq.(81) and taking into account that 
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• fc-point correlations with k - odd vanish, 

• the /c-point correlations factorize into a product of the correlation functions of cor- 
responding subsequences, as soon as the distance between any two points exceeds 
unity, 

we find that Z n {k) are polynomial functions of t of order [n/2], where [x] is the floor 
function, of the following form 

Zn(k)=[^-) £(-1)^, (82) 

in which expression W^ n denote the weights of corresponding configurations. These weights 
Wi in can be obtained via straightforward combinatorial calculations and are given explicitly 
by: 

w,, n = e (- - 21 + 1 ) f s r: j ) ! , n (c-t . <83) 

\ Z^i=i m j J rni\m2 ] .m 3 \ . . . mi'. -M- 
for / > 0, where the summation extends over all positive integers mj obeying the equation: 
mi + 2m 2 + 3m 3 + . . . + lm\ = I, and 

Wi=0,n = 1. (84) 

Note that in Eq.(83) the symbol g), in a standard fashion, denotes the binomial coefficient 
for a > b and equals zero otherwise. 

Now, it is expedient to introduce the generating function of the form: 

oo 

Z(k,z) = J2Zn(k)z n . (85) 

n=2 

Note that summation over n is performed for n > 2, since Af = for n — 1 and n = 0. 
Multiplying both sides of Eq.(82) by z n and performing summation, we get, (relegating 



22 



the intermediate steps to Appendix B), 

4 2 



Z(k,z) = 



+ 



'l + & k fz l + e lk 
4 



z + 



(1 + e ik fz 



(86) 



Therefore, Eq.(86) relates the characteristic function of the sum of peaks and throughs in 
random permutations, (or, equivalently, of the number of turns of the PGRW trajectories), 
to the /c-point correlations in the rise-and-descent sequences, Eq.(74). Using next the 
expression in Eq.(73), we find, eventually, the following explicit result for the characteristic 
function Z(k, z): 



Z(k,z) 



{l + e ik f z 



l_ e i*\V2 



1 + e 



ik 



coth((l-e-) 1/2 0-l 



1 + e 



ik 



-z. (87) 



Expanding next Z(k, z) in Taylor series in powers of variable k, we get 

- 2 2z 2 ik z 2 (60 + I5z + 6z 2 - z 3 ) k 2 



Z(k,z) 



+ 



(4 



l-z ' 3(1-^) 2 90 (l-^) 3 2 

which implies, in particular, that the average and averaged squared number of turns of the 
PGRW trajectory are given by 



2(n- 1) 



(89) 



and 



(5n 2 - 7n + 2) 



(n-3) 



(n 2 -7n + 12) 



|2 fi( n _ i) + ^_^e(n - 3) + ^ ^(n-4). (90) 



Consider finally the form of Z(k, z) in the limit of small k and z — > 1 .In this limit, we 
find from Eq.(87) that in the leading order Z(k, z) obeys: 

1 



Z(k,z) 



2i , 4 ,2 



(91) 



Inverting Eq.(91) with respect to /c and z, we thus obtain that in the asymptotic limit n — > 
oo, the distribution function V(Af,n) of the number of "turns" of the PGRW trajectory 
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converges to a normal distribution of the form: 



3/5\i/2 



P(A/»~- 

4 \nn 



cxp 



V 



45( W- -n 

16n 



/ 



(92) 



This is not, of course, a counter-intuitive result in view of the polynomial three-term 
recursion in Eq.(76) obeyed by (closely related to V(Af,n)) peak numbers P(n + l,m) [23]. 



7 Diffusion limit 

Consider finally a continuous space and time version of the PGRW in the diffusion limit. 
To do this, it is expedient to define first some sort of an "evolution" equation for Vi(Y). 
We proceed as follows. Define first the polynomial: 

V«\x,Y)=Y,Qy 1 (x), (93) 

Y l= Y 

where the polynomial Qy l (x) has been determined in Eq. (64) and the sum extends over all 
/-steps trajectories Y\ starting at zero and ending at the fixed point Y . Note next that one 
has 

V l (Y) = (vW(x,Y)) {ti} (94) 

Now, for the polynomials V^(x,Y) one obtains, by counting all possible trajectories Y\ 
starting from zero and ending at the fixed point Y, the following "evolution" equation: 

V {l+1) (x, Y)=i r V {l \x, Y-1) + I r V {l \x, Y + 1). (95) 

Taking next advantage of the established equivalence between the processes Y x and X^, 
we can rewrite the last equation, upon averaging it over the distribution of variables {xi}, 



as 



- (1 i^r^ Y - + w^) nY + 1] ' (96) 

which represents the desired evolution equation for Vi(Y) in discrete space and time. 
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We hasten to remark that Eq.(96) can be thought of as a direct consequence of the 
celebrated relation between the Eulerian numbers: 



We emphasize, however, that despite the fact such a relation is evident for the distribution 
Vi(X® = X) of the end-point of the PGRW trajectories, Eq.(ll), its generalization to 
the distribution of the intermediate steps Vi(xj n ^ = X) as well as to the distribution of 
the process Y h Vi(Y = X), is a completely unevident a priori result, which enlightens the 
process under consideration. 

Given the evolution Eq.(96), we turn next to the diffusion limit. Introducing space 
y = aY and time t = rl variables, where a and r define characteristic space and time 
scales, we turn to the limit a, r — > 0, supposing that the ratio a 2 /r remains fixed and 
determines the diffusion coefficient D = a 2 /2r. In this limit, Eq.(96) becomes 



Note that the resulting continuous space and time equation is of the Fokker-Planck type; 
it has a constant diffusion coefficient and a negative drift term which, similarly to the 
Ornstein-Uhlenbeck process, grows linearly with y, but the amplitude of the drift decays 
in proportion to the first inverse power of time, which signifies that the process Y\ eventu- 
ally delocalizes. Note also that the form of Eq.(98) ensures the conservation of the total 
probability. 

The Green's function solution of Eq.(98), remarkably, is a normal distribution 



which is consistent with the large-n limit of the discrete process derived in Eq.(18). 
Now, the Langevin equation corresponding to Eq.(98) reads 




(97) 




(98) 




(99) 



dy 

dt 



y 
t 



+ C(t) 



(100) 
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with (((t)((t')) = 2D 5(t — t'). Its solution can be readily obtained and has the following 
form: 

y(t) = ]J\>((t')dt' (101) 

Note that Eqs.(98) or (100) model dynamics of an overdamped particle in a one-dimensional 
continuum subject to a force — y/t and a white noise. Note also that Eq.(lOl) implies (see 
Ref. [24]) that y{t), being a linear functional of a white noise, can be represented as a 
rescaled Brownian motion of the form: 

y{t)= (\'{t')dt' (102) 



where {C(t)('(t')) = 2D 5(t — t')/3. Consequently, y(t) can be also thought of as the 
solution of the customary Langevin equation with zero external force, 

f t =C(t). (103) 

We finally remark that Eq.(103) could have been derived directly from the definition of 
the discrete process Y h by taking the proper continuum limit of Eq.(55). The correlator 
(C(t)C(t')) must be taken in this case as the limit of the continuation C^ 2 \m) of C^ 2 \m) 
to m < 0: 

(C'(f)C'(0> = lim n C (2) M = I™ (s mfi - \{8 m ^ + S m:1 )) = 2D 5(t - 0/3 (104) 
The Gaussian nature of the process y(t) is transparent within this formalism. 



8 Conclusions 

In conclusion, we have studied here a simple model of a non-Markovian random walk, 
evolving in discrete time on a one-dimensional lattice of integers, in which the moves 
of the walker to the right or to the left are prescribed by the rise- and- descent sequence 
characterizing each random permutation 7r of [n + 1] = {1, 2, 3, . . . , n + 1}. We have 
determined exactly the probability V n (X) of finding the end-point X n of the walker's 
trajectory at site X. Furthermore, we have shown that in the long-time limit V n (X) 
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converges to a normal distribution in which an effective diffusion coefficient D = 1/6 is three 
times smaller that the diffusion coefficient (D — 1/2) of the conventional one-dimensional 
Polya random walk. This implies that correlations in the generator of the walk - random 
permutations, which arise because only a finite amount of numbers is being shuffled and 
neither of these number may be equal to each other, are marginally important. Indeed, 
we have shown that two- and -four-point correlations in the sequence of rises depend 
only on the relative distance between them and extend effectively to nearest-neighbors 
only. Next, at a closer look on the intrinsic features of the PGRW trajectories, we have 
formulated an auxiliary stochastic Markovian process. We have shown that, despite the 
fact that this process is Markovian, while the random walk generated by permutations is 
not a Markovian process, the distribution of the auxiliary stochastic process appears to 
be identic to the distribution of the intermediate states of the walker's trajectories. This 
enabled us to obtain the probability measure of different excursions of the permutations- 
generated random walk, determine general /c-point correlation functions and to evaluate 
the asymptotic form of the probability distribution of the number of turns in an n-step 
PGRW trajectory. Finally, we have discussed, in the diffusion limit, the continuous space 
and time version of such a walk. 
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10 Appendix A 



Here we briefly outline the derivation of the integral representation of V n (X) in Eq.(15). 
Using the following expansion of the polylogarithm function [16]: 



q= — oo 



q + 



7Y 



n+2 



we have that V n (k) in Eq.(29) can be cast into the form 



7T 



q=— oo 



q + 



7Y 



n+2 



Noticing next that 



E 



q= — oo 



1 



q + 



7Y 



n+2 



(n+l)\ 

(-l)" +1 7T n+2 d n+1 



*n + l ( 1 - - )+ (-l)^n + l (~ 



7T 



7T 



cot(fc), 



(105) 



(106) 



(107) 



(n+1)! dk n+1 

where ^ n {x) denotes the polygamma function [16], we arrive eventually at the representa- 
tion in Eq.(15). 
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11 Appendix B 



Here we present the details of the derivation of the generating function in Eq.(86). Multi- 
plying both sides of Eq.(82) by z n and performing summation, we get 

(1 + e ik ) 



1=0 v 7 n=2« V 



1 + e ife 



1 + e ifc ^ \1 + e ifc / I ^ mi!m 2 !m 3 ! . . .mj! v ; 

oo 



Z = 



X 



^ fn - 21 + 1\ (z(l + e ik ) 



l-e 2 



-V z 2 ) V v J y TT(c (2j) ) 



i=0 x \ .7=1 

+ 1 \ (1 + e ik )Y 2 



x 



Note now that ( E ( + i 1 m .) = for p < Ej = i m j ~~ 1- Hence, for I > 0, (when, evidently, 
Y^j=i m i — 1)> we have 



y,/ p + 1 \ z(l + e ik ) Y ~ / p+1 ^(l + e* fc ) 



z_(l + e*)\ Si = imj_1 / _ g (1 + e- fe ) ^ "^^^ 



(109) 



On the other hand, for / = (when V ■ = 0), 



g^^y.^^y'. (110) 

Making use next of Eqs.(109) and (110), as well as of the identity 

E - z E , n V ' , II *7 = r^i — i - ( m ) 

f— ' ^ mi!m 2 !m 3 ! . . . m,! J i _ \ x r j 
we arrive eventually at the result in Eq.(86). 
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